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Abstract. Recent models of the envelopes of seven mas- 
sive protostars are used to analyze observations of Hjj~ 
infrared absorption and H 13 CO + submillimeter emission 
lines toward these stars, and to constrain the cosmic-ray 
ionization rate £cr- The H 13 CO + gives best-fit values of 
Ccr= (2-6 ± 1.8) x 10 -17 s , in good agreement with dif- 
fuse cloud models and with recent Voyager/Pioneer data 
but factors of up to 7 lower than found from the H^ data. 
No relation of £cr with luminosity or total column den- 
sity is found, so that local (X-ray) ionization and shielding 
against cosmic rays appear unimportant for these sources. 
The difference between the and H 13 CO + results and 
the correlation of A r (H^") with heliocentric distance sug- 
gest that intervening clouds contribute significantly to the 
Hg" absorptions in the more distant regions. The most 
likely absorbers are low-density (< 10 4 cm -3 ) clouds with 
most carbon in neutral form or in CO. 
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1. Introduction 

The ionization fraction of molecular clouds is an important 
parameter for their dynamics through its control over the 
influence of any magnetic field. The ionization also has 
a major effect on the chemistry of molecular clouds be- 
cause ion-neutral reactions are generally much faster than 
neutral- neutral reactions. In dense regions shielded from 
direct ultraviolet irradiation, the ionization is dominated 
by cosmic rays. However, the rate of this process £cr has 
not yet been constrained directly. The current best esti- 
mate comes from chemical models to reproduce the ob- 
served abundances of OH and HD in diffuse interstellar 
clouds (Hartquist et al. 1978; van Dishoeck & Black 1986; 
Federman et al. 1996), notably those toward Perseus OB2. 
These models indicate that Ccr= 1CP 16 — 10~ 17 s _1 per 
H atom, but with a factor of 10 uncertainty because of un- 
certainties in temperature, radiation field and the effects 



of shocks. In addition, it is unknown if £cr varies with 
location in the Galaxy, since the diffuse cloud results are 
limited to the solar neighbourhood. 

Each cosmic-ray ionization of H2 yields one H^ 
molecule, so that H^ has a constant concentration 
which depends only on £cr and the abundances of 
its main destroyers: CO, O and electrons. Hence, 
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van der Tak et al. (2000) of the envelopes of these stars, 
and use the derived values to model observations of 
HCO + , the abundance of which is also proportional to 
Ccr- 

2. Models 

McCall et al. (1999) present observations of rovibra- 
tional lines of H^" in absorption against seven lumi- 
nous (10 4 — 10 5 Lq) young stars, which are still em- 
bedded in envelopes of ~ 100 M Q of dust and molec- 
ular gas. These same sources have been studied by 

13 CO infrared absorption and 



Mitchell et al. (1990) 

by van der Tak et al. (2000)| in submillimeter dust con- 
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tinuum and CS, C 34 S and C 17 line emission. Based on 
these data sets, van der Tak et al. (1999, 2000) modeled 
the temperature and density structure of the envelopes 
using a power law structure n = no(r/ro)~ a . The ra- 
dial dust temperature profile is calculated self-consistently 
from the luminosity and no is determined from submil- 
limeter photometry which probes the dust column density. 
The parameter a is constrained by modeling the relative 
strengths of the CS and C 34 S J = 2 -> 1 through 10^9 
lines with a non-LTE radiative transfer program based on 
the Monte Carlo method. 

The outer radii of the models are twice the half- 
intensity radii of the CS J — 5 — > 4 emission, given in Ta- 
ble 5 of van der Tak et al. (2000)| For the sources studied 
here, the values are (3 — 11) x 10 17 cm, which are accu- 
rate to a factor of 2. However, the dust and gas maps also 
reveal emission extending outside the envelopes of W 3 
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Table 1. Observed and modeled column densities and ortho/para ratios of H^" and CO. 

Source Observed (a) Modeled (6) Inferred Ccr (c) 



JV(H 3 +) o/p N(CO) AT(H+) o/p JV(CO) (1) (2) 

10 13 cm" 2 ratio 10 19 cm" 2 10 13 cm" 2 ratio 10 19 cm" 2 KT 17 s _1 



GL 2136 


38 ±4 




2.2 


3.9 


0.45 


3.6 


9.7 


16 


GL 2591 


22 ± 2 


o.8i° 9 ; 2 


1.3 


2.0 


0.46 


2.1 


11 


18 


GL 490 


11 ±6 




0.78 


5.2 


0.19 


8.4 


2.1 


23 


W 33A 


52 ± 13 


0.8i° 1 


2.6 


16.6 


0.26 


7.8 


3.1 


9.3 


NGC 2264 


< 12 




2.2 


3.1 


0.29 


7.2 


< 3.9 


< 13 


W3 IRS5 


< 8.6 




2.6 


4.3 


0.51 


9.4 


< 2.0 


< 7.2 


S 140 


< 4.4 




0.74 


2.3 


0.41 


6.3 


< 1.9 


< 16 



H+ from |McCall et al. 



(1999); C O from 13 CO observations by Mitchell et al. (1990, 1993, 1995). 



Using the physical structure from van der Tak et al. (2000) and assuming £cr = 1 x 10 s 



Case (1): (iV b s (H+) /iV model (H+)) xlO" 



(2) = case (1) x(iV obs (CO) /iV model (CO)). 



IRS5 and GL 490, and maybe W 33A. No "skin" appears 
to surround GL 2591 and GL 2136. The case of S 140 is 
complicated: the dust appears to be heated by multiple 
sources, and its emission is not well fitted by a centrally 
heated model. We do not have a dust map of NGC 2264, 
but this source is part of an extended molecular cloud 
complex, so that extended material may also contribute 
to the absorption. Evidence for extended components 
at lower temperature and/or column density than those 
probed by the dust emission comes from emission in low- 
J lines of CO at > V offsets, and self-absorptions on their 
central line profiles. These features are present in the data 
for all the sources discussed in this paper. Before consid- 
ering foreground contributions in § |^, we concentrate on 
the dense molecular envelopes. 

3. Results 

Given the temperature and density profiles, we calcu- 
late the H^" concentration at each position in the en- 
velopes. Considering only cosmic rays as producers of Hg 
and reactions with CO and O as destroyers, the concen- 
tration of Hjj" is given by n(R^)— CcR/[a ; (CO)xfcco + 
x(0)xko]- In this expression, a;(CO) and x(0) are the 
abundances of CO and O relative to H 2 , and kco and 
fco the rate coefficients for their respective reactions with 
H3 , taken from Millar et al. (1997). We neglect any de- 
pendence of fcco on temperature since the dipole moment 
of CO is small. The models use an abundance of CO of 
2 x 10~ 4 at temperatures above 20 K, and zero below due 
to freeze-out on dust grains. This abundance behaviour 
is consistent with observations of C 17 emission lines by 



van der Tak et al. (2000). The abundance of O is assumed 
to be 1.5 x 10~ 4 based on the models of Lcc et al. (1996) 



and the temperature in our models does not drop below 
14 K, where O would freeze out. The ortho/para (o/p) or 
(J, K) — (1, 0)/(l, 1) ratio of Hg changes with radius since 
the ground state of ortho-H + lies 32.86 K above that of 



para-Hg ( Dinclli et al. 1997 ), and reactive collisions with 



H2 tie the o/p ratio to the kinetic temperature. Figure [l] 
illustrates the results for the case of GL 2136. 

Integration of the concentrations of ortho- and para- 
Hjj" over radius yields total column densities ./V(H;j~) and 
mean o/p ratios, which are compared with the data in 
Table |l|. The calculated o/p ratios are consistent with the 
data within the observational errors, but the model values 
are systematically lower than those observed. For CO, the 
models, which were constrained by emission data, typi- 
cally overproduce absorption measurements of iV(CO) by 
factors of 3, probably due to deviations from spherical ge- 
ometry on small scales, consistent with several other trac- 
ers (van der Tak et al. 2000). Since the model values of 
A^Hj) may be less affected because Hjj~ is more evenly 
distributed than CO (Fig. 1), Table |lj presents the values 
of Ccr both before (case 1) and after (case 2) scaling the 
model down by the ratio of observed to modeled iV(CO). 
The uncertainty in the model is a factor of two due to 
the uncertain radii of the envelopes. The estimates of £cr 
in case (2) are considerably larger than those in previ- 
ous work (§ [l]), which together with the low o/p ratios 
indicates that there may be an additional component of 
warm along the line of sight. We will estimate the con- 
tribution to A^(H^) by the dense envelopes by modeling 
emission lines of H 13 CO+ which have critical densities of 
~ 10 6 cm' 3 , and hence cannot arise in the foreground. 

4. Comparison with HCO+ 

In the dense envelopes, the main destruction route of 
is the reaction with CO into HCO+. The concentration of 
HCO+ is given by n(HCO+)=x(CO)n(H+)fc C o/ [{x(e)k e 
+ .T(H20)fcH 2 oL with k e the rate coefficient for dissocia- 
tive recombination of HCO+. The electron fraction x(e) 
has been calculated at each point in the envelopes with 
a small chemical network (cf. de Boisanger et al. 1996) 
based on the UMIST reaction rates (Millar et al. 1997). 
The main difference with the analysis of de Boisanger et 
al. (1996) is the use of a detailed physical structure to in- 
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Table 2. Observed and modeled fluxes 
(K km s- 1 ) of the H 13 CO+ J = 3 -> 2 and 4 - 



/ T mh dV 
3 lines. 



Source 


Observed 


Modeled 


CCR 


iV(H+) 




3-2 


4-3 


3-2 


4-3 


(a) 


(b) 


GL 2136 


4.1 


4.0 


4.5 


2.9 


3.3 


8 


GL 2591 


5.5 


4.4 


5.5 


3.2 


5.6 


7 


GL 490 


2.4 


1.9 


3.0 


2.0 


0.64 


0.3 


W 33A 


8.9 


5.3 


9.8 


5.3 


1.3 


7 


NGC 2264 


4.0 


2.1 


3.3 


2.4 


0.61 


0.6 


W3 IRS5 


9.3 


9.2 


8.2 


6.0 


2.8 


3 


S 140 


9.2 


8.6 


8.5 


6.6 


3.7 


1.1 


(a): Best fit 


value of Ccr to H 


13 CO+ 


data, in 


lO" 17 B- 


-l 



(6): Predicted from H 13 CO + for case (2), in 10 13 cm" 



terpret the high-excitation lines. We assume that O2 and 
H2O have negligible (<10~ 6 ) abundances in the bulk of 
the envelopes, but that at T > 100 K, ^(H^O) jumps 
to 5 x 10~ 5 due to grain mantle evaporation. We neglect 
metals such as Mg, Fe and S as contributors to x(e) and 
large molecules such as polycyclic aromatic hydrocarbons 
as sinks of x(e); using the low metal abundances inferred 
from dark cloud chemistry models would increase x(e) by 
a factor of 2-3 (Lee et al. 1996). The values of (cr derived 
above give x(e) ~ 10~ 7 at the outer radii and ~ 10~ 9 at 
the inner radii, as illustrated in Fig. [l]. The precipitous 
drop of HCO + at 100 K, caused by reactions with evapo- 
rated water, occurs at too small radii to affect our results. 

In the comparison with data, we use the 60 x less abun- 
dant isotope H 13 CO + to avoid optical depth effects. The 
maximum optical depth in the lines is ~ 1 in our models. 
Table I lists the calculated fluxes of the H 13 CO+ J=3^2 
and 4^3 lines in 18" and 14" beams. Observations are 
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Fig. 1. Temperature and density structure of GL 2136, 
and calculated concentrations of fTj", e and HCO + , both 
with (full line) and without (dotted line) destruction by 
H 2 included. 



from |van der Tak et al. 



(1999)| for GL 2591 and from 



le Boisangcr et al. (1996)| for NGC 2264 and W 3 IRS5. 
The data for W 33A, GL 490, S 140 and GL 2136 were 
obtained with the James Clerk Maxwell Telescope in the 



way described in van der Tak et al. (2000) . Using £cr de- 



rived from H^" , the models overproduce HCO + by factors 
of 2 - 7. Adjusting the models to the H 13 CO + data yields 
refined estimates for (cr (Table ||) which pertain strictly 
to the dense molecular gas, unaffected by any intervening 
clouds along the line of sight. The data for the various 
sources span the range of Ccr= (2-6 ± 1.8) x 10~ 17 s _1 , 
in good agreement with the diffuse cloud estimates (§ Q), 
and also consistent with recent data from the Voyager and 
Pioneer spacec raft at distances up to 60 AU from the Sun 
( IWebber 1998b - 

Figure |2| shows that the source-to-source variation in 
Ccr is not related to Galactic structure through differences 
in cosmic-ray flux, nor to shielding against cosmic rays at 
high H2 column densities. The values of £cr are also un- 
related to luminos ity, which implies th at local ionization 
such as by X-rays ( Maloney et al. 1996 ) is unimportant on 
the scales traced by our data. Variations of the cosmic-ray 
density by 50% on scales of a few kpc are in good agree- 
ment with results from 7— ray observations (e.g., Hunter 
et al. 1997). However, why does H^" give systematically 
higher values of Ccr? 



5. Contributions by foreground layers 

Figure [| plots the observed iV (Hj") versus heliocent ric dis 



tance, including all data from McCall et al. (1999) as well 
as the results for the Galactic Center and the diffuse cloud 

The 



(1999) 



in front of Cyg OB2 #12 from peballe et al. 
dense cloud data have a correlation coefficient of 93%, 
suggesting that absorption by intervening clouds plays an 
important role for the more distant sources. This section 
investigates the possible nature of these absorbers. 
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Fig. 2. Derived cosmic-ray ionization rates versus galac- 
tocentic radius (left) and NCH2) in a 15" beam (right). 
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Fig. 3. Observed iV^Hg") versus heliocentric distance. The 
correlation suggests that intervening clouds are important. 



First, the absorptions may occur at the edges of 
the dense cores studied here, where carbon is in 
neutral or ionized form. This "photodissociation re- 
gion" occupies ~3-4 magn itudes of visual extinction 
( Hollcnbach fc Ticlcns 1997 ), corresponding to Nh~8 X 



10 cm . The ionized layer is negligible because the high 
electron fraction (~ 10~ 4 ) limits «-(H^") to 10~ 7 cm -3 . For 



the neutral component, assuming n ~ 10 cm as derived 
specifically for S 140 by fr immermann et al. (1996)], and 



few xlO 13 cm 2 



n(H+) - 10~ 4 cm" 3 , we find N(B^ 
comparable to the dense envelopes. 

Second, the absorbers may consist of cold ( < 20 K) 
molecular gas. For n>10 4 cm~ 3 , the CO will be frozen 



out on the grains. Tielens et al. (1991) observed solid CO 
in absorption toward all our sources and found 7V(CO)~ 
10 17 cm -2 , or TVh ~ 10 21 cm -2 , assuming that most car- 
bon is in solid CO. The implied column lengths are too 
short to be of importance for H;j~, and the low tempera- 
tures are incompatible with the observed o/p ratios. 

Third, the absorptions may arise in clouds with 
n^lO 4 cm -3 , which either surround the power-law en- 
velopes or happen to lie along the line of sight. Such 
"translucent" foregrounds are visible in our data (§ |J), 
and can co ntribute iV(H^) 



models by van Dishoeck & Black (1989). These tenuous 
clouds have long path lengths and may dominate the H^ 
absorption. At low densities, HCO + may form through 
OH + C+ -> H + CO+ followed by H 2 + CO+ -> H 
+ HCO + . However, for our sources, [CII] 158 (im data 
indicate N(C + )/N(GO) <10~ 2 . Translucent clouds are 



generally weak in HCO + emission ( Gredel et al. 1994 ). 



The velocities of the H^" absorptions are consistent 
with those of the submillimeter emission lines of C 17 and 
C 34 S, suggesting that the H^ absorbers are in the vicinity 
of the infrared sources. However, the correlation of iV(Hg ) 
with distance remains after subtracting the dense core 
contribution (Table ||), suggesting a non-local origin. Al- 



together, the data indicate that the contribution of the en- 
velopes to iV(H^" ) varies from ~ 10 13 to ~ 10 14 cm~ 2 , and 
that any additional absorption seen in sources at d > 2 kpc 
occurs in intervening clouds. 

In summary, observations of H^ absorption and 
H 13 CO + emission lines, combined with models of the 
temperature and density structure of the sources, con- 
strain the cosmic-ray ionization rate to Ccr= (2.6 ± 1.8) x 
10~ 17 s _1 , with upper limits that are factors of 3-5 higher. 
Future tests of the results include more sensitive observa- 
tions of H+ toward W 3 IRS5 and NGC 2264, velocity- 
resolved observations to search for H^ absorption at off- 
set velocities from the dense cores, observations of more 
distant sources to test the correlation with distance, and 
observations of HCO + infrared absorption lines to directly 
compare with H;j" and CO infrared absorption. 
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